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ABSTRACT 

Difference imaging provides a new way to discover gravitationally lensed 
quasars because few non-lensed sources will show spatially extended, time vari- 
able flux. We test the method on lens candidates in the Sloan Digital Sky Survey 
(SDSS) Supernova Survey region from the SDSS Quasar Lens Search (SQLS) 
and their surrounding fields. Starting from 20768 sources, including 49 SDSS 
quasars and 36 candidate lenses/lensed images, we find that 21 sources including 
15 SDSS QSOs and 7 candidate lenses/lensed images are non-periodic variable 
sources. We can measure the spatial structure of the variable flux for 18 of these 
sources and identify only one as a non-point source. This source does not display 
the compelling spatial structure of the variable flux of known lensed quasars, so 
we reject it as a lens candidate. None of the lens candidates from the SQLS sur- 
vive our cuts. Given our effective survey area of order 0.71 square degrees, this 
indicates a false positive rate of order one per square degree for the method. The 
fraction of quasars not found to be variable and the false positive rate should 
both fall if we analyze the full, later data releases for the SDSS fields. While 
application of the method to the SDSS is limited by the resolution, depth, and 
sampling of the survey, several future surveys such as Pan-STARRS, LSST, and 
SNAP will avoid these limitations. 
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Introduction 



Gravitational lensing has many applications, from exoplanet searches to large scale 
structure (see the reviews by Kochanek, Schneider, and Wambsganss in the Saas Fe lec- 
tures, 2006). Galax y scale lenses can be used to study the dark matter mass profile of 
the lens galaxy (e.g. Kochanek I Il99ll ; iRusin and Kochanek I l2005t iKoopmans et al. 1 12006 



Jiang and Kochanek 



l2007h . When the background source is a quasar, microlensing by 



individual st ars in the lens galaxy can be used to probe the structure of t he quasar's ac 



cretion disk (Poinde xter. 



regio ns ( Eigenbrod et al. 



ters (IRefsdal 



1964 



Vlorgan. and Kochanek 1120071 ; iMorgan et al. 1120071 ) and broad line 



2007). Galaxy scale lenses can also constrain cosmological parame- 



Oguri Il2007al ). Unfortunately, most applications require large samples of 



gravitational lenses to be competitive with other methods, while fewer than a hundred lensed 
quasars are known. Moreover, these lenses were discovered using different methods with dif- 
ferent biases, which is especially p roblematic if hom ogeneous samples with well-understood 
selection functions are needed 



sec 



Kochanek 1120061 ) 



The known lenses were found by their morphological structure or the presence of higher 
redshift features in the spectrum of a lower redshift galaxy. Morphological surveys ex- 
amine optical or radio images of quasars for evidence that they are lensed. This works 
best for point-like sources like optical quasars, and flat spectrum radio sources. The two 
largest s earches for lensed AGN are the Cosmic Lens All-Sky Survey of flat spectrum radio 
sources (IMvers et al. M2003I : iBrowne et al. 1120031 ) and the SDSS Quasar Lens Search (SQLS; 
Oguri et al. 1120061 ) . The radio samples are limited by the number of sufficiently bright radio 
sources and the difficulties in obtaining source redshifts. The optical quasar samples suf- 
fer from confusion from stars and galaxies and the effects of color changes, either from the 
starlight or dust in the lenses, on sample selection. The spectrosc opic method made se veral 
serendipitous discoveries of lensed quasars, such a s Q2237+0305 (IHuchra et al. 1119851 ) and 



SDSSJ090334.92+502819.2 (Johnston et al. 



Kochanek 1992 



SLACS survey (IBolton et al. 



Miralda-Escude and Lehar 



2004 



2003) . Following theoretical investigations (e.g. 



19921 : iMortlock and Webster ||2000L l200lh . the 



200 



(]) used the spectroscopic method on massive early- 



type galaxies in the SDSS to identify ~50 candidate lensed starforming galaxies, and con- 
firmed 19 of 29 using Hubble Space Telescope images. The spectroscopic method is limited 
by its low yield (about 1 lens candidate per 1000 luminous red galaxy spectra) and biases in 
its mass range (due to the size of the spectroscopic aperture). Moreover, by selecting targets 
based on the properties of the lens galaxy, these lenses are mainly useful for studying the 
lens galaxies rather than for cosmology (see iKochanek 1 120061 ) . 



Kochanek et al. I (120061 ) proposed a new method to find lensed quasars based on dif- 



ference imaging. Difference imaging, also known as image subtraction, is a way to measure 
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the variable intensity in a region of sky (ITomaney and Crotts Ill996l ; lAlard &: Lupton 111998 
Alard 1120001 ) . It has been used to search for a broad range of v ariable sources , such as planets 



trans iting stars (e.g. lHartman et al. 112004 ). microlenses (e.g. lAlcock et al. Ill999l; IWozniak 



2000h . and supernovae (e.g. the Sloan Digital Sky Survey II Supernova Survey, ISako et al 



2007h . In difference imaging, a reference image is made by averaging a set of the best images 
for a field. Then, for each epoch of observation, a difference image is created by subtracting 
the reference image convolved to the PSF and flux scale of that epoch. If a source is variable, 
there will be a flux residual in the difference image corresponding to the variability of the 
source. 

In difference images, lensed quasars are recognizable because they consist of multiple 
variable images that are close together. Since most quasa rs are variable, with 60% having 
o g > 0.05 variability over two years (e.g. ISesar et al. 1120071 ). each of these images is variable. 
The level of variability is then enhanc ed by rmcrolensing of the quasar images by the stars 
in the lens galaxy (see the review by IWambsganss I [2006J ). Thus, lensed quasars look like 
compact clusters of variable sources or extended variable sources, which allows us to easily 
search for lensed quasars in difference images. The number of false positives from pairs of 
variable stars, pairs of quasars (r elated or not), variable star-quasar blends, and supernovae- 
AGN pairs is expected to be low (IKochanek et al. 11200a ). Candidates found this way can be 
confirmed by light curve analysis, since each image will have the same intrinsic variability, 
but with time delays between the images and so me additional uncorrelated variability from 



microlensing by the stars of the lens galaxy (e.g. iPindor II2005I ). Difference imaging can also 
resolve variable lensed sources blended with non-variable sources, such as the lens galaxies. 
This is important because, as we search for fainter lensed quasars, contamination from the 
lens galaxy becomes a steadily greater problem. Tests of image subtraction on the lens 
Q2237+0305 show that the method can find qu asars even when they a re buried by the flux 



of an extraordinarily bright foreground galaxy (IKochanek et al. 1 120061 ) . 



The f orthcoming, l arge scale synoptic surveys, like LSST (ITyson et al. 1120021 ) and Pan- 



STARRS (IKaiser 1120041 ). will be ideal for this method, since they will cover large areas of the 
sky, sample variable sources fre quently, and have deep magnitude limits . In the meantime, 



the SDSS II Supernova Survey (IFrieman et al. 1 12008c : ISako et al. 1 120071 ) is the best survey 



currently available for our goals. Intended to find supernovae for dark energy studies, the 
Supernova Survey repeatedly images SDSS Stripe 82, a 2.5° wide swath along the celestial 
equator, covering 300 square degrees and stre tching across the Sou thern Galactic Cap from 



right ascension 300° to right ascension 60° (IFrieman et al. 1 120081 ) . Ten to twenty public 



epochs were available for a given field whe n we started this project, with th e number rising 



to about forty in the most recent releases (lAdelman- McCarthy et al. 1 120071 ). 
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We applied the method outlined in iKochanek et al. I (120061 ) to the SDSS Supernova 



Survey fields that contained 26 previously identified candidates for gravitational lenses, with 
39 total components, from the initial SQLS candid ate list. These in clude 15 candidates 
with 24 components in the main statistical sample (llnada et al. 1120081 ). and 11 additional 



candidate systems, with 15 images, that were outside the final selection criteria. These 
candidates are listed in Table [TJ We also checked to see if any pair of lens images in a 
candidate lens appeared on the difference images, since the images can be well-separated. The 
lens candidates had already been rejected for other reasons, such as different spectral energy 
distributions for the postulated lens images, or the lack of a lens galaxy. Unfortunately, there 
are no confirmed lensed quasars in the Supernova Survey region. We had two goals. First, 
to use the variability method as an independent check of the SDSS candidates. Second, to 
get a sense of the false positive rate from any other variable sources in the fields. We outline 
the method, our approach to source selection, and the results when applied to targets in §2. 
We discuss the variable sources in the lens candidate fields in §3, including quasars (§3.2) 
and an extended variable source (§3.3). Finally, we conclude in §4 by discussing prospects 
for finding gravitational lenses in future surveys with difference imaging. 



2. Procedure 



2.1. Preparation and Image Subtraction 



For each field around the lens candidates, we downloaded the r-band images from the 
SDSS Data Server, which had epochs t hrough 2005. We chose r-band because the r-band 
SDSS images have great effective depth (jYork et al. 1120001 ). We then created images centered 
on the candidate's RA by combining the image containing the candidate with either the 
previous or next image in the drift scan with an appropriate RA offset. The RA centered 
images had the same dimensions and pixel scale as the original SDSS images - 2048 pixels 
by 1489 pixels with a pixel scale of 0'.'396. The median FWHM seeing in the images was 3.23 
pixels (1'.'3). We masked data within 10 pixels of bad pixels, because these prevent clean 
image subtraction. The bad pixels were usually caused by saturation on the images, so the 
masked areas differed little from one epoch to the next. Since the blank areas could not 
be used for either the reference image or image subtraction, we merged the masks from all 
observations of a particular field into a common mask. We also enlarged the masks on the 
difference images themselves before creating the "absolute value" images (see §2.2), in order 
to minimize edge effects. We lost 15.8% of the area covered by the images to masking. 



After using ISIS (lAlard 1120001 ) to register the images, we created a reference image. We 
used the nine epochs with the best seeing, as measured by PSF size. If fewer than 18 epochs 
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were available, we instead used only images with better than median seeing. The median 
seeing of the reference images was 2.35 pixels, or 0'.'93. Finally, we created the difference 
images. We used ISIS to align the images, determine the convolution kernel needed to 
convert the PSF of the reference image to that for each epoch, and then carry out the image 
subtraction. In some epochs, the difference image did not seem clean, perhaps because of 
poor seeing or thin clouds. To compensate, we eliminated the two epochs for each source 
with fluxes that were furthest from that source's mean flux, and we did not use images in 
which the measured magnitude of the sources in the field was on average more than 1 a from 
the mean. We limited our structural analysis to sources with peak surface brightnesses \i in 
this image with \i < [iback + 0.5 relative to the local background \iback-, where fi and fiback are 
in magnitudes per pixel area. 



2.2. Photometry 

Once we had the difference images, the next step was to select candidate variable sources 
in each field. We used ISIS to extract the uncalibrated light curves of sources in the field 
of each candidate. While we were not concerned with the absolute flux of the sources, only 
their variability and their relative brightnesses, we did set an approximate calibration using 
SDSS stars (point sources) within 3' of each quasar. We used the average offset between 
the ISIS light curves and the SDSS magnitudes to calibrate each field. We selected our lens 
candidates in three steps. 

First, we made a variability cut on the r—a r (magnitude— magnitude variance) plane. 
We excluded sources with more than 30000 counts (r < 17.1), because they could be satu- 
rated. We also excluded sources with less than 2% variability (cr r < 0.022) to account for the 
survey's photometric accuracy. Finally, we fit a line to the r— log a r relation. Those sources 
that had a a r greater than the 4 a upper bound on the line fit were considered variable. Quan- 
titatively, we considered sources to be variable only if they satisfied log 10 ay > 0.362r — 8.56. 
These cuts are shown in Figured) 

Second, quasars have long term, non-periodic variability (LPV) rather than periodic 
variability, so we identified quasars as objects whose light cur ves are significantly better fit 



as parabolas than constants. Following lHartman et al. I (120041 ) . we used the F-test to gauge 



the significance of the improvement of a parabolic fit over a constant light curve. The F-test 
is based on the ratio 

XV(N - 2) 1 1 

where Xp/{N — 4) is the x 2 P er degree of freedom for the parabolic fit, Xc/{N — 2) is the 
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X 2 per degree of freedom of a constant light curve, and N is the number of epochs being fit. 
We calculated F after dropping the two epochs most poorly fit by a constant lightcurve, and 
those epochs with bad images, as described in §2.1. The expected mean of F, if x\ an d x\ 
are independent chi-square distributions and x 2 c nas (N ~ 2) degrees of freedom, is 

N-2 

and the expected dispersion in F is 

2 _ 2(iV-2) 2 (2iV-8) 

(iV-4)3(iV-6) 1 ' 

so we use 

5, = (4) 

to estimate the significance of the improvement from using a parabola. Since xf> an d Xc 
are not independent, with x\ < x\ f° r eacn light curve, fip is not the actual mean of the 
calculated F. It can be shown that Sp < for all sources. A cut of Sf < —0.7 seemed to 
work fairly well when we examined the light curves by eye, although the choice was somewhat 
arbitrary. Raising the threshold would increase the completeness at the price of introducing 
more false positives, while lowering it has the opposite effect. 

Finally, to evaluate the structure of the variable flux, we summed the absolute values 
of the difference images from each epoch. Variable sources appear as positive peaks in 
these "absdiff" images, whatever their light curves. We then used SExtractor to identify 
and characterize the variable sources in these images, focusing on sources that passed our 
variability and F-test cuts. We measured the effective half-light radius R e containing 50% 
of the flux. We also measured the elongation, defined by the ratio of the major axis to the 
minor axis fit to the source. Sources that had large effective radii (R e > 2.5 pixels = 0'.'99) 
or were elongated (e > 1.75) were considered extended. 



3. Variables in the Fields 

Table [2] lists the number of sources that survive each cut and some alternate combina- 
tions of the cuts. When we consider all sources, after excluding bad images (usually due 
to low transparency) and dropping two epochs from each light curve, we could only ana- 
lyze 23 of the 26 fields - the light curves for the SDSSJ2134-0054, SDSSJ2139-0114, and 
SDSSJ2200— 0107 fields had too few epochs for analysis. In these 23 fields, which cover an 
unmasked area of 0.71 square degrees, we can measure light curve statistics, such as a r and 
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Sf, for 15064 objects. We have statistics for 29 candidate images from 21 candidate lens 



system s of the SQLS, as well as 44 spectroscopically identified quasars in ISchneider et al 



(120071 ). Seven candidate lens images and 14 spectroscopic quasars could be analyzed in the 



"absdiff" images, with only one spectroscopic quasar being extended. 

Figure [1] shows the mean magnitude and the standard deviation of each light curve, 
excluding any dropped epochs. Our sources span the range 14 < r < 24. The variance 
shows the usual dependence on magnitude, scaling with the photon noise until saturation 
effects begin to dominate at r < 15. At the typical magnitude of an SDSS quasar, r ~ 19, 
the typical variance is of order 0.01 magnitudes. Variable sources lie above the r—a r curve, 
and we see that these include some, but not all, of the SQLS lens candidates and the other 
quasars in the field. If we apply the cuts described in §2.1, we find 370 variable sources of 
all kinds. 

In Figure [2J we plot the distribution of the sources in the space of Sf and a r . We 
chose Sf < 0.7 as our cut, since it works well in combination with the cut on the overall 
variability. After making those two cuts, we were left with 21 sources, including eight of the 
SQLS lens candidates and 15 of the SDSS quasars. A quick inspection of the light curves 
shows that the F-test does in fact pick out those light curves with a long term trend. In 
Figured we show four light curves for sources that pass the r— a r cut, three of which also pass 
the F-test. The light curve at upper left ( SDSS J2040 12.68-003040.36) i s one of the sources 



that did not pass the F-test. This s ource is neither a ISchneider et al. I (120071 ) spectroscopic 
quasar nor an SQLS lens candidates ( Inada et al. 20081 ). Note the outlier points, of the kind 



that compelled us to drop some epochs as described in §2.2. While the light curve has some 
variability, especially in the earlier epochs, there is no overall trend, and so it fails the F-test. 
The other three light curves pass the F-test and show trends in thei r variability. These three 



variable sources consist of a non-quasar, an SQLS lens candidate (jlnada et al. II2008I ). and 
a spectroscopic quasar which also seems to be extended. These light curves are for sources 
that either clearly fail the F-test (Sf 3> —1) or clearly pass it (Sf < —1) - for sources near 
the cut (Sf ~ —0.7), the light curves and the F-test results are more ambiguous. Our F-test 
cut is generous, in that most of the sources passing it do not clearly pass (Sp < —1), as can 
be seen in Figure |5J 

A few of the sources that do pass the variability cuts and the F-test cannot be analyzed 
on the "absdiff" images. Eighteen of the 21 sources that passed the variability cut and the 
F-test were measured by SExtractor on the sums of the absolute values of the difference 
images. All pass our surface brightness cut of fi ma x < Hback + 0.5, where [i ma x is the peak 
surface brightness per pixel after background subtraction, and \iback is the surface brightness 
of the background per pixel. We chose this surface brightness as a measure of whether the 
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SExtractor detection was real, or a spurious detection of noise. So, we find that 18 sources 
have significant long period variability and measurable spatial structure to their variability, 
out of the 15064 initial sources with light curves. 

As we survey in Table El we can compare these results with those obtained by omitting 
either the r—a r cut or the F-test. In the former case, we are left with 208 sources altogether 
that pass the F-test. Of theses, 85 were detected by SExtractor on the absdiff images, with 
only 54 passing our surface brightness cut on the "absdiff" images so that they could be 
analyzed for spatial structure. This discrepancy, where some SExtractor detections seemed 
to be spurious, motivated our surface brightness cut. If, on the other hand, we include the 
r—a r cut but ignore the F-test, then out of the 370 variable sources, 311 were detected on 
the "absdiff" images by SExtractor. Of these, 158 passed our surface brightness cut and can 
be analyzed for their spatial structure. The combination of the two cuts is very effective at 
reducing the number of candidates relative to the number of quasars known to be present in 
the fields. 



3.1. The SQLS Lens Candidates 



We looked for both extended SQLS sources and for pairs of sources for the SQLS can- 
didates in the difference images. Of the original 26 lens candidates in the SDSS Supernova 
Survey region, we find that some have variable components, but none are extended. As 
Table [2] shows, only eight passed the F-test, and seven could be analyzed on the "absdiff" 
images. None satisfied the criteria given in §2.2. A small separation lensed quasar will ap- 
pear as a cluster of variable sources on a difference image. Thirteen of the SQLS candidates 
are well-separated image pairs, but in all these candidate lenses, at most one image survives 
the F-test, a nd the component that p asses the tests is always the spectroscopically identified 
quasar from ISchneider et al. I (120051 ) . and not the candidate second component selected by 
color criteria. Since none of the SQLS lens candidates show extended variability, and since 
no pair of candidate images have two variable components, none of them are se lected by our 
method either. We conclude, on independent grounds, that llnada et al. I (120081 ) were correct 
to reject them as lens candidates. 



3.2. Quasars in the Candidate Fields 



Using extended variability to find gravitationally lensed quasars assumes that each 
quasar is in fact variable enough to appear on difference images. To address this ques- 
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tion, we also examined the images of the 60 SDSS spectroscopic quasars ([Schneider et al 
20071 ) in the fields of the candidate lenses, including the SQLS lens candidates themselves. 
As summarized in Tabled 49 of the 60 were detected by ISIS; six of the non-detections were 
in masked regions, and one, a SQLS lens candidate, was blended with an object identified by 
Inada et al. I (120081 ) as a galaxy. We were able to perform the least squares fitting necessary 
for the F-test on 44 quasars. Of these 44, 15 passed our usual variability cuts and had 
Sp < —0.7, and we could measure the spatial structure for 14 of these. So, about one third 
(32%) of the quasars for which we have light curve statistics show variability on the difference 
images. Only one of the quasars in the candidate lens fields, SDSSJ213245. 24+000146. 4, was 
extended using the criteria given in §2. It was not one of the lens candidates in the SQLS, 
and we discuss it in §3.3. 



Sesar et al. I (120071 ) examined the Supernova Survey region and found t hat the majority 



of un resolved quasars listed in an earlier release of the SDSS Quasar Catalog (jSchneider et al. 
20051 ) with g < 20.5 were variable. Specifically, over the two years of the survey, over 60% 
of the quasars had variabilities over a g > 0.05, and over 90% of the quasars had variabili 



ties ov er a g > 0.03. When we consider a r alone, we find less variability than ISesar et al 



(120071 ) . Of the forty quasars with statistics and r < 20.5, we find only eight (20%) have 
variabilities greater than a r > 0.05, and twenty-three (58%) have variabilities greater than 
ay > 0.03. Indeed, Figure [1] shows that many quasars, especially the fainter ones, are not 
highly variable. According to our variab ility cut on r and oy, twenty-five quasars (57%) 



are variable, closer to the proportions in ISesar et al. I (120071 ). but still significantly below 
90%. The discrepancy i s partly caused by our use of r-band images, since the variability i s 



greater in bluer bands (ICutri et al. I Il985t IVanden Berk et al. 1 12004 ; IWilhite et al. 1 1200a ) , 
and partly by the shorter time coverage of our data. 



3.3. Candidate Extended Objects 

Figure |H shows the elongations e (ratios of major to minor axes) and effective (half-light) 
radii (R e ) of the sources that are variable, or have Sf < —0.7. Most sources cluster around 
R e ~ 2 pixels; since the pixel scale for SDSS is O'.'396/pixel, this is consistent with a seeing 
of about Of! 8. Most also have elongations of e < 1.5. There are some outliers, that are either 
elongated (e > 1.75) or big (R e > 2.5). For comparison, we measured the seeing of each 
epoch and compare the R e distribution to the HWHM = FWHM/2 distribution of the data. 
For a Gaussian PSF, the HWHM encloses 50% of the light. As can be seen in Figure IU 
most of the images have measured seeing less than our R e cut, with a median seeing of 
1.6 pixels (0'/6) and 6.5% having HWHM > 2.5 pixels. Only one source passed all the cuts, 
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SDSSJ213245.25+000146.5. 



Unlike most variable sources, which only marginally pass the F-test, SDSS J213245. 25+000146. 5 
clearly passes it with Sp = —1.117. The difference images themselves (Figure [5] bottom 
far- right), however, show no clear structure in the variability. Furthermore, this source only 
marginally passed our R e cut, with R e = 2.8 pixels, and it is not elongated. This suggests that 
i t is unlikely to act ually be a lensed quasar. SDSSJ213245. 25+000146. 5 is an X-ray source 
dVoges et al. 2000 ) and one of the quasars in the SDSS Quasar Catalog ( Schneider et al. 



20071 ). Since quasars are in fact variable, it is not surprising that it should pass our r— a. 
cut. However, the quasar is at a redshift of only z = 0.2342, while typical le nsed quasars are 



expec ted to be more distant, with essentially none at such low redshifts (e.g. iMitchell et al 



20051 ). In addition to showing the quasar, the ref erence image (as in F i gure ED seems to show 



the host galaxy of SDSSJ213245.25+000146.5. IVanden Berk et al. I (120061 ) considered the 



source to be resolved and were able to fit galactic components to the spectrum, finding that 
2% of the flux was from the apparent host galaxy. If SDSSJ213245. 25+000146. 5 is a lensed 
quasar, one might expect there to be two se ts of spectral lines a t different redshifts as in the 



spectroscopic method of finding lenses (e.g. iBolton et al. II2006I ). 



4. Conclusion and Future Prospects 



We searched for sources with the spatially extended variability characteristic of gravita- 
tionally lensed quasars by applying difference imaging to the Sloan Supernova Survey fields of 
the SQLS lens candidates and their surroundings. We found one source, SDSSJ213245. 25+000146. 5, 
that passed basic criteria for variability, non-periodic variability, and which appeared ex- 
tended on its "absdiff" image. However, it did not show the com pelling structure seen in 
difference images of known lensed quasars (IKochanek et al. 1 120061 ) - it only barely passed 
the effective radius criterion, suggesting it is not likely to be a true lensed quasar. Further- 
more, it was not among the SQLS lens candidates in the SPSS Supernova Survey region. 



Our criteria successfully rejected the SQLS candidates, as had llnada et al. I (120081 ) using 



othe r criteria. If used with other lens search techniques, like color or morphological crite- 
ria (jOguri et al. 1120061 ). difference imaging could substantially reduce the number of false 
positives, although it would only work for quasars that are variable during the observation 
period. Only a third of the quasars detected by ISIS passed our quasar variability tests 
and could be analyzed for spatial structure in the difference images. More generally, as a 
variability survey for lensed quasars covering 0.71 square degrees, we successfully identified 
one third of the SDSS quasars in the field and had no false positive detections of lenses. 



In many ways, our experiment was limited by the data. First, there are no known lensed 
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quasars in Stripe 82, so we could evaluate no examples of success or set any limits on the 
role of false negatives. Furthermore, when we carried out this experiment, the Supernova 
Survey had few public (< 25) epochs available. Since some epochs will be dropped because 
of bad seeing, or because image subtraction failed, there are some fields where there were 
too few enough epochs to extract light curve statistics. Similarly, the released epochs did 
not cover a large enough time interval. In some cases (e.g. SDSSJ213245. 25+000146. 5 in 
Figure [3]), all the non- rejected epochs spanned only a few months. The performance of our 
approach would improve markedly with survey data spanning several years and with many 
more epochs. Additionally, the relatively poor resolution of the SDSS survey is not ideal 
for id entifying; qua s ar len ses. This is also true for color and morphological techniques, such 



as m llnada et al. I (120081 ). where candidate lenses with image separations of less than l'/O 
were rejected. With better resolution and more epochs, false positives in the method would 
decrease. If we examined every SDSS quasar in the supernova region, we would expect to 
find of order 1000 were variable with ~ 50 false positives if we assume no improvement in 
the false positive rate with the addition of more epochs. In reality, we should have a marked 
reduction in the rate because we could focus on epochs with better seeing and add in the 
variability information from the ugiz bands as well. 

While these shortcomings will be partly solved by the full SDSS Supernov a Survey data 



set, t he real future for th e method lies with ground-based surveys like LSST (ITyson et al 



20021 ) and Pan- STARRS (IKai 



scr 



2004) and 



Aldering 



20021 ) 



proposed surveys from space like the SuperNova 
Pan-STARRS and LSST would repeatedly 



Acceleration Probe (SNAP, 
survey very large areas (10 3 -10 4 rather than 10 2 square degrees) with im proved resolution 



jp'.h—l ". rather than > l'/O) and depth (~ 24 mag rather than ~ 22 mag). iKochanek et al. 
( 120061 ) estimated that LSST can discover ~ 10 3 lensed quasars with V< 23. The space based 
surveys would cover far less area (15 square degree for SNAP) but with vastly improved 
resolution (~ 0'/05) and depth ( ^ 28 mag). Sin c e few lenses have separations this small 
(l'/5 is the expected median, e.g. iMitchell et al. I (120051 )). any lens identified by variability 
is easily confirmed by its morphology. More important for the space missions is that their 
great depth allows searches for other lensed, time variable sources such as the supernovae 
themselves. 
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Table 1. SQLS LENS CANDIDATES IN THE SDSS SUPERNOVA SURVEY REGION 
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Component 1 is always a spectroscopically identified quasar given in ISchneider et al 



(120051 ). If the candidate images are less than 2'.'5 apart, the components are blended, and no 
position is given for Component 2. If they are more than 2'.'5 apart, Component 2 is selected 



-16- 



by color criteria and the position is given. 

b Positions from i nitial SQLS c a ndida te list, but these candidates were outside the final 
selection criteria for llnada et al. I (120081 ). 



'Positions listed in llnada et al. I (120081 ). 
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Table 2. NUMBER OF SOURCES SATISFYING EACH CUT 



Criteria 


ISIS sources 
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a As in the initial SQLS candidate list (llnada et al. II2008I ). 
b As given by lSchneider et al. I (120071 ). 



c Determined by r-oy cut listed in §2.1, for sources with statistics 
d Determined by Sp < — 0.7 
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c Determined by /i < ^back + 0.5 

f Determined by e > 1.75 or R e > 2.5 pixel 
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Table 3. CANDIDATE EXTENDED VARIABLE SOURCES 



Object Field m a m Sf e -R e a 



Pass r—a r Variability Cut 
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a Given in units of pixels. The image scale is 0'/396 pixel 1 . 



-20- 



10 



0.1 



□ 01 - 



~i 1 r~ 









0.001 



0.0001 



_i i i_ 



_L 



_L 



ISIS source 
O 3QLS lens candidate 
• 3DSS spectroscopic quasar 
Elongated or large 



source (from Table 3) 



_L 



_L 



14 



16 



IB 20 
r (mag) 



22 



24 



Fig. 1. — The standard deviation of the light curves versus average magnitudes. Our vari- 
ability cut is marked by the black line segments - sources above the boundary are considered 
to be variable. 
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Fig. 2. — The F-test significance Sf, as described in §2.3, plotted against the standard devi- 
ation of the magnitudes. Our cut of Sp < —0.7, indicating long term (LPV-like) behavior, 
is drawn on the plot as a solid line. Most of the extended candidates only pass it marginally, 
with Sf > —1.0 (dashed line). 
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Fig. 3. — Light curves of variable sources passing the r—a r cut with different Sp values. 
Epochs that were cut are marked in red, with crosses meaning the epoch was on average 
more than 1 a away for all sources, and rings meaning the epoch was one of the two furthest 
from the mean for that particular source. The sources o n top are not known q uasars, while 



the two on the bottom are SDSS spectroscopic quasars ([Schneider et al. 11200 



_qu 

3). 
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Fig. 4. — Elongation e and effective radius R e for sources with peak surface brightness 
A* < /Ufeacfc + 0.5 on the absolute value images, which we considered to be real detections by 
SExtractor. Our cut of e > 1.75 or R e > 2.5 pixels is also plotted; those sources above or to 
the right of the cut were considered extended. On the bottom, we show a histogram of the 
seeing (FWHM/2) of the original images. 



-24- 




Fig. 5. — Various sources on the reference images (top) and the "absdiff" images (bot- 
tom). The images are centered on the objects, and all have the same scale. At 
far left is SDSSJ205201.63+000051.6, a quasar with low variability (a m = 0.019 and 
Sf = —0.20), which is not detected by SExtractor in the "absdiff" image. The two 
middle sources (SDSSJ021247.22+003433.5 at mid-left, SDSSJ021328.83+003030.8 at mid- 
right) pass the F-test and are extended, but do not pass the r—o r cut. At far-right is 
SDSSJ213245. 25+000146. 5, our one candidate extended variable source. 



